Exact coherent states of a harmonically confined Tonks-Girardeau gas 
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Using a scaling transformation we exactly determine the dynamics of an harmonically confined 
Tonks-Girardeau gas under arbitrary time variations of the trap frequency. We show how during 
a one-dimensional expansion a "dynamical fermionization" occurs as the momentum distribution 
rapidly approaches an ideal Fermi gas distribution, and that under a sudden change of the trap 
frequency the gas undergoes undamped breathing oscillations displaying alternating bosonic and 
fermionic character in momentum space. The absence of damping in the oscillations is a peculiarity 
of the truly Tonks regime. 



One of the most challenging and interesting problems 
in quantum dynamics involves understanding a tempo- 
ral behavior of strongly correlated many-body systems 
beyond the linear-response regime or the adiabatic ap- 
proximation. Aside from fundamental interest, this issue 
is of primary importance for current experiments with ul- 
tracold atomic gases and their potential applications to 
quantum information, where control and manipulation 
of entangled states of many-particle quantum systems is 
required. 

The one-dimensional (ID) gas of impenetrable bosons 
(Tonks-Girardeau gas) corresponds to the limit of in- 
finitely strong repulsive interactions in the Lieb-Liniger 
(LL) model of ID bosons interacting through a contact 
pair potential It was shown in Ref. [2| that the 
Tonks-Girardeau limit is applicable for describing the 
low-density regime of bosonic atomic gases in a quasi- 
1D geometry. The physical properties of ID bosons in 
this limit can be investigated in detail since their wave- 
function is known explicitly in terms of the one of a non- 
interacting fermions in the same external potential 0] . In 
fact, the density profile, the thermodynamic properties, 
the collective excitation spectrum and the density cor- 
relation functions coincide with those of an ideal Fermi 
gas, leading to interesting manifestations of fermioniza- 
tion of a Bose gas, such as broadening of the density 
profiles [4, an increase of the frequency of collective ex- 
citations ya] and a dramatic reduction of the three-body 
recombination rate 0. However the one-body density 
matrix, and consequently the momentum distribution, 
differs considerably from that of a Fermi gas, due to the 
phase correlations stemming from the bosonic statistics 
of the Tonks gas. For the homogeneous case the mo- 
mentum distribution n(p) at the origin has a peak 
and for a harmonically trapped gas the population 
of the lowest single-particle state scales as VN with N 
being the particle number 0, Q ; this shows that due to 
the strong interactions the bosons do not form a Bosc- 
Einstein condensate. In both cases at large momenta 
p > hn, with n being the density at the center, the mo- 
mentum distribution shows characteristic slowly decay- 
ing tails n(p) ~ p~ 4 [Io|. 

Experiments on cold atomic gases under optical con- 



finement in a quasi-one dimensional geometry are now 
starting to explore the strongly interacting regime, 
demonstrated by the examination of its correlation prop- 
erties llj and of the frequency of collective modes 
More recently there has been a significant progress to- 
wards the Tonks limit, in which the momentum distribu- 
tion ^3 and the thermodynamical properties have 
been measured. Experiments addressing the dynamics of 
the Tonks gas seem also in view. Some aspects of the dy- 
namical evolution of a Tonks gas have already been theo- 
retically studied. These include the formation of solitons 
in a ring geometry [lij . the splitting and recombination 
of a Tonks beam across an obstacle and ID expan- 
sion I n the last case, using numerical calculations 
on a lattice, the momentum distribution was found to 
approach that of an ideal Fermi gas . 

The aim of this work is to study the dynamical evo- 
lution of a harmonically trapped Tonks gas induced by 
arbitrary time variations of the trap frequency. We show 
that this evolution can be described exactly using time- 
dependent coherent states (see e.g. 0), in close analogy 
of the dynamics of Bose-Einstein condensates [2(| • These 
have been investigated also for the case of a strongly cor- 
related Bose gas interacting with an inverse-square pair 
potential |2lj . Here we explore in particular how the 
bosonic or fermionic properties of the Tonks gas mani- 
fest themselves in the dynamics of the coherence. For this 
purpose we choose the momentum distribution as the ob- 
servable. We first determine explicitly its time evolution 
in terms of the initial-time configuration with the aid of a 
scaling transformation. We then use this result to study 
two important examples: (i) the ID expansion of the gas, 
where we explain why the gas develops a Fermi shape of 
the momentum distribution, and (ii) the large-amplitude 
breathing modes, where we find the absence of damp- 
ing and a rich dynamical evolution in momentum space, 
displaying alternating bosonic and fermionic character. 

Scaling transformation. We consider N impenetrable 
bosons with mass m in a ID geometry at zero tempera- 
ture and subjected to a harmonic potential 

V ext (x,t) =mu; 2 (t)x 2 /2 (1) 

with u>(t) arbitrary time-dependent trapping frequency 
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and ui(t < 0) = Wo- We proceed to derive an 
exact analytic expression for the evolution of the 
Tonks gas wavefunction, the one-body density ma- 
trix and the momentum distribution at all times. 



We start by employing the time-dependent Bose-Fermi 
mapping |l5| . which allows to write the many-body 
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wavefunction $y(a;i, x%, . . . , x^; t) for the Tonks gas 
in terms of the wavefunction <&p(xi, X2, ■ ■ ■ , XNl~t) — 
(1/ViV!) detj fc=1 (f>j(xk,t) of a noninteracting Fermi gas 
experiencing the same external potential. The single- 
particle orbitals <j>j(x,t) satisfy the time-dependent 
Schrodinger equation 
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2m dx 2 



(t>j(x,t) + Vea : t(x,t)(f>j(x,t). (2) 



The Tonks gas wavefunction is then constructed 
by applying to the unit antisymmetric function 
A{xi, . . .,x N ) = TLi<j <k < N sgn(xj - x k ), 



$t(xi, . . .,x N ;t) = A(xi 



,x N )^ F (x 1 ,...,x N ;t), 



(3) 

and is thus properly symmetrized. 

For the case of a time-dependent potential Q the 
introduction of a scaling transformation for both the 
space and time coordinates provides an exact solution 
for Eq. ©, which reads [l9ll20| 



<f>j(x,t) 
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) exp 
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(4) 



In Eq. (@J above the scaling factor b(t) obeys the second- 
order differential equation 



b + L0 2 {t)b = UJ 2 /b 3 



(5) 



with initial conditions 6(0) = 1 and 6(0) = 0, the rescaled 
time parameter is determined by r(t) = j Q dt 1 /b 2 {t'), and 
<j>j(x, 0) are the well-known wavefunctions of the ID har- 
monic oscillator with frequency ojq and eigenvalue Ej ex- 
pressed in terms of the Hermite polynomials. We remark 
that Eq. is the unique time-dependent solution of the 
linear Schroedinger equation ©. 

Substituting the one-particle states (3$ into Eq. Q 
leads to the final result for the time evolution of the Tonks 
gas wavefunction in terms of its initial-time expression 

$ T ( Xl , ..,x N ;t) = b- N / 2 <t> T ( Xl /b, ..,x N /b; 0) 



xex p(^Ei|] ex p(- j E^ 



with Iq — A^/ft/mwo- Here we used that (i) the scaling 
(0} takes place irrespectively of the one-particle quantum 
number, and (ii) the unit antisymmetric operator A is 
invariant under the scaling transformation. 

Equation © allows an immediate derivation of an 
analytic expression for the one-body density matrix, 



given by its first-quantized expression gi(x,y;t) = 
N Jdx2--dxN&*{x, X2, •■, xn] t)Q(y, X2, ■■, xn] t), in the 
scaling form 



gi(x,y;t) = ^gi |;0^ 



exp 



i b x 2 — y 2 



buj 



2l 2 



(7) 



This result shows first of all that during the generic time 
evolution described by Eq. (JIJ the absolute value of the 
one-body density matrix preserves its power-law behavior 
at large distances [9|, |23| ■ For the case of an expansion 
this was noticed in Ref. |18| from numerical simulations 
on a lattice. Secondly, Eq. J7J) yields the exact evolution 
of the Tonks gas density profile p(x;t) = qi (x, x\ t) = 
p(x/b;0)/b, thus generalizing the result of Ref.|lj| for 
any time-evolution of the trap frequency. We remark that 
Eq. J7J also holds at any finite temperature, as follows 
from the statistical Bose-Fermi mapping theorem j2^ . 

The dynamical phase in Eq. J2J enters the expres- 
sion for the time evolution of the momentum distribution 
n(p,t) = J dxdy e lp ^ x ~ v ^ h gi(x,y;t), which upon change 
of integration variables x/b, y/b — * x, y reads 



(8) 



n{p,t) = b / dxdy gi(x,y;0) 



x exp 



-ib 



b x 2 — y 2 p(x — y) 



w 
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Below we describe two cases where the dynamical phase 
acquired by the one-body density matrix strongly influ- 
ences the time behavior of the momentum distribution. 

Expansion. A ID expansion could be achieved in an 
experiment by turning off only the longitudinal confine- 
ment. To describe this case we set cu(t < 0) = uoq and 
oj(t) — for t > 0, and the solution of Eq. © for the 
scaling parameter b(t) is given by b{t) = yl + uj^t 2 . As 
this parameter becomes increasingly large with time, we 
are able to determine analytically the behavior of the mo- 
mentum distribution (JHJ for long times by the method of 
stationary phase. The points for which the phase is sta- 
tionary are given by x* = y* — loqpIq / '(6ft) , and hence we 
find that the momentum distribution is asymptotically 
determined solely by the diagonal part of the equilibrium 
one-body density matrix, i.e., the particle density profile 
p(x; 0) which is identical to the one of an ideal Fermi gas. 
For the case of harmonic confinement the latter is propor- 
tional to the equilibrium momentum distribution uf(j>) 
of the ideal Fermi gas |24| thus allowing to write the final 
expression for the Tonks gas momentum distribution as 



n(p,t) = |w /6| n F (u> p/b). 



(9) 



This result can be also understood as follows. Since the 
initial momentum distribution of the atoms is very nar- 
row, the final momentum distribution is determined by 
the hydrodynamic velocity field acquired during expan- 
sion. Such a velocity field is obtained from the dynamical 
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FIG. 1: Momentum distribution of an expanding Tonks gas 
with N = 7 at different times (in units of 1/ljo) as indicated on 
the panels, from numerical solution (solid lines), asymptotic 
fermionic limit (dotted lines) and Thomas-Fermi approxima- 
tion (dashed lines). The units are indicated on the axis labels. 
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FIG. 2: Momentum distribution of an oscillating Tonks gas 
with N = 9 and ujo/uji = 10 at different times (in units of 
T — 7r/wi) indicated on the panels, from numerical solution 
(solid lines) and Thomas- Fermi approximation (dashed lines). 
The units are indicated on the axis labels. 



phase in the Eq.JJJ) and is linear in the position. Hence 
its distribution has the same shape as the particle density 
profile of a Fermi gas. For N 3> 1 Eq. @ is well described 
by the Thomas-Fermi approximation, which neglects the 
quantum shell oscillations of order 1/N |2f| . 

n(p,t) ~ 2{ll/%)\u Q /b\ yjp* - (wqp/6) 2 . (10) 

with Pp = hV2N I lo being the Fermi momentum. 

We estimate the characteristic time, tp, for such 
"dynamical fermionization" by considering that in the 
Thomas-Fermi regime the one-body density matrix de- 
pends on coordinates only through the dimensionless ra- 
tios x/R and y/R. By rewriting Eq. @ in such rescaled 
coordinates and in terms of p/Pf we see that the large 
parameter governing the dynamic phase in Eq. JHJ is 
Nbb/uio, and thus we find tp ~ 1/Nluq ~ h/Ep, Ep 
being the Fermi energy. 

Equation (PJ) provides an accurate description of the 
momentum distribution for long times and for small mo- 
menta p < Pp. We proceed now to derive the second 
exact result regarding the large momentum behavior of 
n(p, t) at all times. From the scaling solution © it fol- 
lows that during its time evolution the many-body wave- 
function displays the same type of cusp singularity as in 
its equilibrium configuration. Since no additional singu- 
larity is present in the dynamical phase, this cusp de- 
termines alone the large momentum behavior of the mo- 
mentum distribution, as in the case of the equilibrium 
solution ^(J- Hence, we find a power-law decay of the 
momentum distribution at large p, with an additional 



time-dependent suppression factor which originates from 
the dilatation of the interparticle distances during the 
expansion, 

n(p,t)^p- i b-\ (11) 

This complements the result (jHJ showing that at time t > 
1 /loo the fermionization of the momentum distribution is 
complete and the large-p tails are negligible. 

We have tested the above predictions and explored 
the expansion at early times by numerically evaluating 
the momentum distribution Eq. (|SJ with a Fast Fourier 
Transform method using the explicit expression of the 
initial time one-body density matrix in terms of a de- 
terminant of Hankel type |9|. The results, reported in 
Fig. ^ show that a broad Fermi distribution rapidly de- 
velops during expansion, followed by small adjustments 
of the quantum shell oscillations, while the p~ A tails be- 
come less and less important at long expansion times. 

Oscillations. As a second example, we consider now 
the case of an abrupt change of the trap frequency which 
induces large-amplitude "breathing" oscillations in the 
gas, i.e. we set uj{t < 0) = u>o and uj(t) = uii for t > 0, 
with lui < ujq- The solution of Eq. (0) is 

b(t) = I + - ul) sin 2 (cj^/ljI (12) 

which describes periodic oscillations between 1 and loq /u>\ 
with period T = ■k 

The dynamical evolution of the cloud in coordinate 
space is described, according to Eq. (7J, by the self- 
similar breathing of the density profile with a time law 
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given in Eq. (|12|) . Notice that the solution i|12|) implies 
that the oscillation is undamped. The time evolution 
in momentum space, given by Eq. ©, displays a richer 
structure, and in particular an oscillating behavior be- 
tween a bosonic-like and fermionic-like momentum distri- 
bution which is illustrated in Fig. [21 The main features 
of the dynamics may be understood by the following an- 
alytical considerations. When the condition Nbb/u>o > 1 
holds, the stationary phase method can be employed, 
yielding that a time-dependent Fermi-like structure de- 
velops as in Eq. ©, and that the large- wavevector tails 
are suppressed as in Eq. In the regime JV > 1 

and <x>o 3> wi this result is valid for most values of t, and 
closely resembles the dynamics of an ideal Fermi gas. The 
latter is given by the exact expression [i^ 
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n(p,t) = B(t)n F (B{t)p). 



(13) 



a small time interval woAi = N^oJq/ (uSq — uf) around 
the turning point t = T/2 the fermionic description does 
not hold. There, using Eq. (|5J) we find that the atomic 
cloud recovers the initial bosonic momentum distribu- 
tion rescaled by a factor b max = uq/oji, according to 
n(p,T/2) = b 

max 

In conclusion, we have shown that the dynamical phase 
acquired by the many-body wavefunction during the evo- 
lution is responsible for the "dynamical fermionization" 
found in the ID expansion and that during a large- 
amplitude breathing mode the cloud displays oscillations 
in momentum space between a Fermi-like and a Bose-like 
structure. 

In this work we have considered only the Tonks- 
Girardeau limit of the Lieb-Liniger model. In this case 
the scaling transformation is exact and predicts a co- 
herent, undamped motion of the cloud, which is a very 
peculiar feature of the Tonks regime. We suggest that 
the absence of damping in the breathing modes may be 
used to characterize the truly Tonks regime: for finite val- 
ues of the coupling strength we expect instead a damped 
motion of the oscillations. Our predictions directly ap- 
ply to the experiments on ultracold atomic gases in tight 
optical traps, where the time evolution of the momen- 
tum distribution should be experimentally accessible by 
allowing a 3D expansion of the cloud 131. The possi- 
bility of manipulating coherently the Tonks gas and the 
absence of damping may make this system interesting for 
the applications to quantum information. 
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with B = b I ' \l 1 + b 2 b 2 / 'io 2 being related to the scaling \J. 
of the kinetic energy (B(t) = 1 for the expansion). In 
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